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Abstract 

We report on QCD radiative corrections to the process pp — * Zbb in the approx- 
imation in which the b quark is considered massless. The implementation of this 
process in the general purpose Monte Carlo program MCFM is discussed in some de- 
tail. These results are used to investigate backgrounds to Higgs boson production in 
the ZH channel. We investigate the Higgs mass range (100 GeV < mn < 130 GeV) 
for the Tevatron running at ^/s = 2 TeV. 



1 Introduction 



In this paper we report on the calculation of the strong radiative corrections to the process 

p + p^Z + b + b. (1) 

These results are obtained from a Monte Carlo program which allows us to obtain pre- 
dictions for any infra-red safe variable. Since the decays of the Z are included we can 
perform cuts on the transverse momenta and rapidities of the final state leptons as well 
as on the properties of the jets present in the event. This discussion is similar to the 
calculation 

p + p + b + b , (2) 

presented in ref . |]] , although there are some differences. Ref . |J can be seen as a first step 
towards the calculation of a vector boson plus 2 jets at 0(a^), where the underlying par- 
tonic process contains only quarks in the initial state. When considering the production 
of Zbb we must also consider initial states containing two gluons, a further stepping-stone 
towards a general V+ 2 jet program. With this in mind, we shall present our method 
in some detail, with a general outline in Section 2 and further details presented in the 
Appendix. 

Much effort has been devoted to the study of Higgs production at the Tevatron running 
at yfs = 2 TeV @. These studies indicate that, given enough luminosity, a light Higgs 
boson can be discovered at the Tevatron using the associated production channels WH 
and ZH. In Section 3 we perform an analysis in the ZH channel that incorporates as 
many of the backgrounds as possible at next-to-leading order. Whilst we use no detector 
simulation and do not attempt to include backgrounds due to particle misidentification, 
the results presented here can provide a normalization for more experimentally realistic 
studies. This is of importance since more detailed studies are often performed using 
shower Monte Carlo programs which can give misleading results for well separated jets. 



2 



2 Calculational overview 



In this section we will outline the implementation of the matrix elements in our Monte 
Carlo program. We separate the discussion according to the three types of contribution: 
the lowest order (Born) processes, the virtual (loop) corrections and the real corrections 
associated with additional soft or collinear radiation. 

2.1 Born processes 

The Born processes which we are considering are 

q + q — > Z + b + b 

g + g^Z + b + b (3) 

which are illustrated in Fig. |I[ For the qq case, the set of diagrams represented in (a) is 
the same as that previously implemented for Wbb production in ref. The diagrams 
in (b) are related to those in (a) by crossing. Diagrams (b) are absent in the Wbb case 
since a W does not directly couple to a bb pair. For the same reason, the gg process (c) 
is also not present in the case of Wbb production and these diagrams introduce a new set 
of matrix elements. In all of the matrix elements we will use the massless approximation 
for the 6-quarks. We expect this to be a good description for bb pairs produced with 
a large invariant mass, m 6 g. This expectation is borne out by Fig. |2| where lowest order 
predictions with and without the 6-quark mass are compared. As expected the corrections 
are of order 4mj/m?r. Note that Fig. [2] may overestimate the mass effect since a fixed 
mass, rather than a running mass which is smaller at high scale, is used. The basic jet 
cuts used in Fig. || to define the bb jets are, 
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Figure 1: Representative diagrams contributing to production at lowest order. A 
complete gauge-invariant set of diagrams is included in the actual calculation. 

where as usual 



The total cross-section and the contribution from ^-initiated processes are shown sepa- 
rately, binned according to the bb invariant mass. If we examine the total cross-sections 
(integrated over a large range, 10 GeV < m b i < 160 GeV), we find that the gluon- 
gluon contribution is 36% of the quark-antiquark result. However, when we are later 
interested in examining this process as a background to a ZH(-^> bb) signal, we will be 
interested in the cross-sections over a limited range of m b i close to a Higgs mass mu, 
where 100 GeV < mjj < 130 GeV. We can see from the figure that in this region the 
gg contribution is a much more significant fraction of the total, roughly 60% of the qq 




(5) 



Fig. |2] also illlustrates the relative importance of the two sub-processes in Eq. (|3]). 
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Figure 2: Lowest order predictions showing the validity of the zero mass approximation 



cross-section in a window around m# = 110 GeV. The reason for the importance of the 
gg process in this region can be understood from the contributing diagrams. By exami- 
nation of Fig. |l|(c) we see that, in contrast to the other diagrams in Fig. [j], m b i may be 
large whilst keeping all propagators on-shell. 



2.2 Real corrections 

The real matrix elements of Nagy and Trocsanyi (NT) given in ref. P] are implemented 
using a subtraction procedure which follows closely the treatment of Ref. ||. To 
illustrate our method, we outline the subtraction terms that are needed for one particular 
piece, namely the leading colour contribution to the gg process. The actual numerical 
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calculation that we have implemented includes the full matrix elements, not just this 
leading colour piece. 

We shall label the momenta of the partons as follows, 

g(-Pi) + g(-P2) -> Z(p u ) + 6(p B ) + 6(pe) + g(pr), (6) 

where all momenta are considered outgoing. In leading colour, it is sufficient to consider 
one ordering of the gluons along the fermion line, and to obtain the remaining terms by 
permutation of the gluons [f|. The real matrix elements are given by, 



l-Mrcail 2 = 8eVWV 3 x( J2 |A(5* 1*2*7* 6*)|%0(i 

i (6 perms) 




= ]T \M 5 (l g ,2 g ,7 g )\ 2 + 0(C F N), (7) 

(6 perms) 

where we have used the (NT) notation for the tree-level amplitude A (see Eqs. (NT:A43- 
A49)Q) and for simplicity we have suppressed all couplings of the Z boson and the sum 
over quark and lepton helicities. As usual Cp = 4/3 and N = 3 is the number of colours. 
The neglected terms are subleading in N. Eq. ([!]) is very similar to the leading-order 
term, 

IA^loI 2 = 8eV^iV 2 x( £ |m(5*l*2*6*) 

\(2 perms) 

= Y |A^4(1 3 ,2 9 )| 2 + 0(C F ), (8) 

(2 perms) 

which we will use to construct the subtraction terms. 
2.2.1 Dipole enumeration 

The matrix elements in Eq. (0) become singular when gluon 7 becomes soft and/or 
collinear with the quarks or either of the other gluons. There are a number of meth- 
ods for cancelling these singularities at next-to-leading order and we choose to use the 
1 We call the reader's attention to the erratum of July, 2000 which corrects some equations in ref. ||. 
The archive version has been updated. 
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dipole subtraction method of Catani and Seymour 0. In this approach, the dipoles pro- 
vide a convenient and efficient way of enumerating the singularities and cancelling them 
in a local fashion, including spin correlations. 

A dipole consists of 3 partons, the first two of which combine to form the emitter 
parton, ai, and the third is the spectator b. The construction of the dipoles starts [|J 
from the eikonal factors present in the soft limit which are rewritten as, 

2p a p b _ 2p a p b r 1 



+ (a <-> b) . (9) 

PaPi PbPi PaPi + PbPi L PaPi J 

The first term on the right hand side corresponds to the soft limit of the dipole with 
emitter ai and spectator b. The full dipole is given by the extension of Eq. @ to include 
collinear emission. The dipole must be labelled by the parton type of both emitter and 
emitted parton, as well as whether the partons lie in the initial or final state. For instance, 
the notation T>\^ refers to a dipole where 17 is the initial-state gluonic emitter splitting 
into gluon 1 and parton 7, with respect to spectator parton 2. 

The dipole structure of the subtraction terms can be understood from the soft gluon 
limit of the matrix elements (|7|). In the soft limit of a single colour ordering (permu- 
tation), gluon 7 can have singularities with neighbouring partons only. Each of the six 
permutations therefore gives rise to two dipoles, leading to twelve terms in all. Hence we 
find, 



l-Mldipoie = N [D^ g + + V^ g + V%; qq + Vi\% g + V%f qq \ x \M*{l g , 2,)|'+(1 <-> 2) 

(10) 

where the dipole functions are listed explicitly in Appendix [A]. 



2.3 Virtual corrections 

The virtual corrections to the processes in Eq. (|3|) are taken from the work of BDK, (Bern, 
Dixon and Kosower For other work on these parton subprocesses, see also ref. [^J. In 
this paper the ultraviolet, infra-red and collinear singularities are controlled by continuing 
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the dimension of space-time to d = 4 — 2e using the four- dimensional helicity scheme. The 
structure of the Monte Carlo program is such that virtual terms are combined with the 
integrals of the dipole subtraction pieces defined above. When added together, e-poles 
cancel and a finite result is obtained. 

To demonstrate this explicitly, we will examine the leading colour contribution to the 
1-loop diagrams. From BDK the leading iV-contribution from the loop diagrams is (cf. 
BDK, Eq. (2.12)), 



Re (2M LO MU oop ) = SeV^iV 3 (g) (4tt) 



x 



4 ree *(5 9 , 1 9 , 2 g , 6,-)A 6;1 (5„ l g , 2 g , 6 ? ) + O 



Extracting only the pole pieces from the leading term (see BDK, Eqs. 
together with the renormalization in Eq. (6.5)) yields, 



+ (l<-2). (11) 
(8.5) and (8.7), 
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+ t +(l<-2), 



where we have adopted the usual definition s^- = 2pi ■ pj. We are now in a position 
to make a comparison with the dipole terms that we presented in the previous section, 
Eq. flTop. Inserting the appropriate integrals (the functions y end from Appendix |A]) we 
find the counterterm from the real contribution to be, 



\M\ 



2 

counter 



X 



+ 



+ 



gig 
^\ 
£l5 

y, 
v. 



+ 



+ 



3e 
3 

2l 
3 

2l 



11 - 

£15 



2n 



N 

-e 

+ 



7T~ 



6e 




2n f \ 
N J 


7T 2 

3 


+ 3 




6e 


("- 


2n f \ 
N J 


3 


+ 3 


+ (!«-► 2 



8 



M Z ,T Z 


91.187,2.49 GeV 


a(M z ) 


1/128.89 


M w , T w 


80.41,2.06 GeV 


Gf 


1.16639 x 10~ 5 



Table 1: Basic input parameters. 



When we add these two contributions together, we see that the e-poles cancel, leaving a 
finite term that is a combination of logarithms and constants multiplying the lowest order 
squared matrix elements, 

Re (2Ai LO A<U op ) POlC + l-^lcounter = 

3 - In (£*U * In -2). (13) 

4 \s 1B J 4 \s 2 J 2 2 J 1 ; 1 ; 

This result exemplifies the cancellation that occurs throughout our calculation, both at 
sub-leading colour and in the other qq initiated sub-process (the details of which we have 
omitted here for brevity). In general, when the poles present in || are added to the 
integrated dipoles given in Appendix [A], the result is extra finite terms proportional to 
the lowest order matrix elements. 



(^)|A<(l.,2,)|" X 



3 Results 

3.1 Basic features 

The standard model is specified by three gauge couplings; the top quark mass also leads 
to large corrections to tree graph results. For our purposes we will replace these four 
parameters by Mw, M z , a(M z ) and Gf, the values of which are given in Table 0. The 
coupling constants e, gw and sin 2 9w are derived from these input parameters according 
to the definitions below, 

e 2 = Ana(Mz) 
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g 2 w = 4V2G F M^ 

sm 2 e w = 4-- ( 14 ) 

9w 

When denned in this way, these are effective parameters which include the leading effects 
of top quark loops [|l(| . In addition, we use the MRS98 parton distribution set with 
a s (M z ) = 0.1175. 

We apply our results to the phenomenological study of the Zbb background to ZH 
production. We consider the channel Z — > vv which has a branching ratio of about ~ 20%; 
the signature in this mode is a 66-pair and missing transverse energy, which we denote by 
Ej>. For convenience we impose the following cuts on our Monte Carlo simulation, and 
these will apply to all the results given in this section. We first require the observation 
of a b and a b jet well separated from each other and from the direction in the transverse 
plane of the missing energy. The cuts we impose (in addition to the basic jet cuts of 
Eq. (D) are, 



(15) 



where <pEy T ,b is the azimuthal angle between the 6-jet direction and the direction of the 
missing Et- We also reject events which have additional jets in the observed region 

bjetl < 2 - 5 > 

|p[ et | > 15GeV. (16) 

We first examine the distribution of the cross-section as a function of the bb mass. Our 
results are presented at leading and next-to-leading order for a renormalization and fac- 
torization scale fi = 100 GeV and also for a much smaller choice fi = 20 GeV in Figs. [3] 
and This smaller scale is chosen because, as can be seen in Fig. ^, for this scale 
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choice the distributions at LO and NLO are comparable for 100 < m bb < 160 GeV. We 
see that there is both a significant enhancement in the total cross-section and a change 
in the shape of the distribution when including the radiative corrections. Moreover, for 
the lower choice of scale, the total cross-section in the range 20 < m bb < 160 GeV is even 
larger and there is a steepening of the distribution towards the peak at low m bb . 

To further investigate the issue of scale dependence in this process, we consider an 
integral of the cross-section over a restricted range of m bb . We will consider such integrals 
in our analysis of the Higgs search, where we shall integrate around the mass of the 
putative Higgs boson. In particular, we will integrate over 

m H - a/2A < m bb <m H + V2A (17) 

where A is the mass resolution of the 66-pair. We assume that A = 0.1 x tuh and perform 
a gaussian smearing with this resolution. The scale dependence of the cross-section is 
shown in Fig. [5[ For low values of the renormalization and factorization scale ft the NLO 
cross-section peaks; it is equal to the lowest order result at fi « 15 GeV. 

To conclude our discussion of the general features of the Zbb calculation at NLO, 
we assess the relative importance of the different sub-processes. We have already seen 
that at leading order the gg process provides a significant contribution to the cross- 
section, particularly for higher values of m bb . For this reason, in Fig. |6] we show the 
cross-section binned by m bb and separated into the contribution of the gg component of 
parton luminosity and the contribution of the rest (beyond leading order we also have qg 
initial states). 

This figure shows that for the interesting range, m hb ~ 100 GeV, the gg contribution 
is even more important, relative to the other sub-processes, than at leading order. In 
fact, if we divide the next-to-leading order curve by the leading order one to obtain the 
'^-factor' we find that for m hb > 60 GeV the ratio K q q/ qg « 1.7 whilst K gg « 2.5. 
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m u T t 


175, 1.4 GeV 




sin 2 9 W 


0.22853 


v ud 


0.97500 


Hig 


gs mass (br) 


100 GeV, (0.8119) 


v us 


0.22220 


Hig 


gs mass (br) 


110 GeV, (0.7697) 


v cd 


0.22220 


Hig 


gs mass (br) 


120 GeV, (0.6778) 


v cs 


0.97500 


Hig 


gs mass (br) 


130 GeV, (0.5254) 



Table 2: Additional parameters used for the signal and background calculations for the 
Higgs search. 



3.2 Higgs search 

Previous studies 0, [12], [T3| have shown that the channel ZH where the Z decays in 
neutrinos (observed in the detector as missing transverse energy) and the Higgs boson 
decays to bb can give a significant signal for the standard model Higgs boson with a mass 
below 130 GeV. 

The largest background to this signal is the Zbb process that we have discussed above. 
After this QCD background, the largest other background at 2 TeV is from the diboson 
process Z°(^ vv)Z{-^ bb). Significantly smaller backgrounds come from the processes 
W ± *{-^ t{— > bW + )b) and q't(—> bW + )b, where in each case the W boson decays leptoni- 
cally, but the charged lepton is not observed. In addition contributions to the signal are 
present from pp — > WH where the jy-decay lepton is missed; correspondingly there are 
contributions to the background from pp — > Wbb and pp — *> WZ. 

In our analysis, we will calculate the significances using our parton level Monte Carlo 
MCFM in which the signal and largest backgrounds are calculated beyond the leading 
order. In order to calculate the signal we require the Higgs branching ratio into bb pairs, 
which is a strong function of the Higgs mass. The values of the branching ratio for the 
four values of the Higgs mass that we will study, as well as additional parameters used to 
calculate the other backgrounds, are shown in Table Q. In order that our normalizations 
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81 


g't(-> 6W+(^ e + i/))6 x 4 


35.6 



Table 3: Values of the total cross-sections for normalization purposes 



are clear, in Table |3] we show the total cross-section for each process, together with the 
results when appropriate branching ratios and lepton flavour sums are included. The 
values shown for the Wbb and Zbb processes only have the basic cuts, Eq. @; for tt, WH 
and ZH the rates shown are the total cross sections with no cuts. 

For each choice of the Higgs boson mass, the signal and backgrounds are integrated 
over an m bb masss range as in Eq. (P7|). For the backgrounds we also reject events with 
additional observed leptons, 

M < 2, 

\pf\ > 19GeV. (18) 

Our results at 2 TeV are given in Table f|, where we have used a double 6-tagging 
efficiency e bb = 0.45. When we compare these results to those of the SUSY-Higgs study [0 
we immediately see a difference in some of the channels. In particular, by including the 
radiative corrections to the Zbb process we have increased this background considerably. 
Furthermore, the contribution to the signal from the WH process and all the backgrounds 
involving a W are significantly lower than in || . The reason for this discrepancy is clear. 
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0.6 


0.6 


0.6 


0.6 


q't{-^ bW + )b 




0.2 


0.2 


0.2 


0.2 


Total B 




39.0 


31.9 


23.6 


19.3 


S/B 


0.17 


0.15 


0.14 


0.11 


s/Vb 


1.04 


0.87 


0.66 


0.46 



Table 4: Signal, backgrounds and significance for the Z-channel at \/s = 2 TeV 
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In our Monte Carlo, all these channels involve an unobserved lepton, which must therefore 
have a very low or be at high rapidity. With a detector simulation, such as is used in 
the SUSY-Higgs study, these unobserved leptons may be central and thus produce large 
contributions to the cross-sections, whilst being unobserved by virtue of misidentification. 

To illustrate the effect of these differences on the significance of the different analyses, 
in Fig. [7] we compare our MCFM results with two approaches from ref. || . The apparent 
similarity of the significances plotted in Fig. [J] is the result of two effects. We believe 
that the significance of ref. || is too high because, the Zbb and ZZ backgrounds are 



underestimated in the light of the work presented here and in Ref. ||15|| . On the other 
hand, our significances are certainly too high, because due to the lack of a full detector 
simulation, the backgounds from processes with mis-identified leptons are underestimated. 



4 Conclusions 

We have presented the first results on the strong radiative corrections to the Zbb pro- 
cess. The results indicate large radiative corrections which can significantly change the 
estimates of the backgrounds to the process pp — > ZH at the Tevatron. We find that the 
significance S/ yB could be less than the result in the report of the SUSY-Higgs working 
group. The results described in this paper do not represent a full analysis of the poten- 
tial to see the Higgs signature in the ZH channel at the Tevatron. The most significant 
shortcoming is the lack of a full detector simulation. 

The primary change is the increase in the size of the Z&fe-background because of 
the inclusion of NLO effects. Our calculation underscores the importance of having an 
experimental determination of the Zbb background, either by relating it to observed Zbb 
events at lower m h i or by relating it to Z+two jet events. 
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A Dipole formulae 



This appendix details the exact dipole functions that we have used in the subtraction 
of the real integral, as well as the integrals of these dipoles that must be added to the 
loop contributions. We note that these functions differ slightly from those given in || 
since they carry no colour factors. In addition, these are the dipoles appropriate for 
the four- dimensional helicity scheme, which is the scheme of choice for the calculation of 
complicated virtual corrections. There are four basic types of dipole, corresponding to 
each of the emitter and spectator partons being in either the initial or final state. In all 
cases we label the partons as ik emitting parton k with respect to the spectator j. 



A.l Initial-initial 

We consider the case of an emitter i and a spectator j, both in the initial state, radiating 
parton k into the final state. We first define the dimensionless variable x to be, 

x = 1 — . (19) 



The dipole functions are then given by, 



via = wl(J—i- x 

xpi ■ Pk V 1 — x 

qq {%) 
XPi-pk 

V%i = (1 - 2x(l - *)) 

xpi-p k 

g 2 H 2e ( 2x 



vZj g = ^—(-^- + 2x(l-x)) (20) 
xpi ■ Pk \ 1 — x J 

These are the dipole functions if we neglect spin correlations between partons % and k 
and are thus not sufficient for the cases T> iitgq and T* ii gg . If the lowest-order amplitude 
is represented by M.^ (-M*), where \i [y] is the polarization index of gluon ik, then the 
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required subtraction is actually, 
where the additional pieces are, 



(21) 
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V 



ik,j 


g 2 ^ 


1 — X 


2pi ■ pj 


ii,gq 


XPi ■ Pk 


X Pi 


■ Pk Pj ■ Pk ' 
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ii, 99 


ii, gqi 








fj, Pi 
= Pk 

Pi 


■Pk fj, 

Pj- 

■Pi 





(22) 



When we perform the integration, it will be over the azimuthal average of these terms. 
The azimuthal average over the tensor q^q v can only depend on the vectors Pi and pj and 
since q ■ pi = it is equal to 



1 Pi- Pi 



(23) 



where A = —q 2 /2. The final result is, 



2 



(24) 



With Vi^gq = V ii:qg = 0. 

To integrate these dipoles we decompose the n-particle phase space in d = 4 — 2e 
dimensions, 

<ty (n W 6 -Pl...Pn)= II \-^-Th 5+ (Pl)} (2*)* ft -f>), (25) 

i=l,n y Z7T ) i=l 

into two pieces, 

rl ~ 

d(/) {n) (p i ,p j ->..., K-i,Pk) = dx d^ (n-1) (pi,Pj ->..., fc„_i) [dp k {pi,Pj, x)} , (26) 

V 

The first term corresponds to an (n — l)-particle phase-space and the second is the one- 
particle dipole sub-space that we must integrate over. Under this decomposition, the 
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momenta of the particles are transformed. The final state momenta all undergo a Lorentz 
transformation (the details, which are given in ref. || erratum, are unimportant here), 
The initial momenta are modified with — > pi = xpi and pj = pj is unchanged. The 
emitted parton phase space is given by, 



[dp k (pi,pj,x)] 



d d Pk 



(2tt)' 



- i 5 + (pi)Q(x)Q(l-x) 



(27) 



Eq. (|27f) can be rewritten as 

[dpk(Pi,Pj,x)} -- 



{2 PiPj y e dn d _ 3 dv dx e(a;(1 _ x)) e(u) @ , 1 



X 



1671-* 



:i -x) 



(27T)!- 



2c 



1 — X 



-2c 



1 - 



.i-*v i-,/J ' (28) 

where x is defined in Eq. (|i~9"D, v = PiPk/PiPj and dQd-3 is an element of solid angle in the 
directions perpendicular to the plane defined by Pi and pj. 

dn d „ 3 (4tt) 



(2.)— r(l-e) =Cr - (29) 

Isolating the common factor g 2 /pvPk from the dipole terms and re- writing all the momenta 
in terms of the transformed ones, we can further express this as, 
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:i-x) 



-2c 



dv 



x 



x 



(30) 



This is now in a form where we can integrate the dipoles in Eq. (pup. We perform the 
integrals and then expand in powers of e, discarding terms of 0(e). In this way we find, 

/^,_\ / ™,,2 

[dpk(Pi,Pj,x)] x{Vfi J ) 



XfX 



2ttJ \2pi ■ pj 



x cr 



1 



Pqq{x) +S(l-X) 



1 



+ 



7T _ 



+ 



41n(l - x) 



e - "• • ' \e* 2e 6 J ' (1 - x) , 

where we have expressed the result in terms of the splitting function 



-2(1 + a;) ln(l-x) 



p qq {x) 



X) 



1 -x + -6(1 -x). 



(31) 



(32) 



Before proceeding further we must perform mass factorization to get into the MS scheme. 
In d = 4 — 2e dimensions, the splitting functions contain additional terms 



Pij(x,e) = Pij(x) + e plf (x) 



(33) 



We have the following expressions for the Pij 

2 



Pqqi 



X 



1 



X) 



3 1 
- 1 - x + -5(1 - x) + e( - (1 - x) + - 5(1 



x , 



Pgq{x) 

pL x ) 



1 + (1 



X 



e x 



X 



x 2 + (1- x) 2 - e 2x(l - x) 



l-x) 



- -A + 2x(l-x) + 5(1 



X 



11 - 2n f /N\ 



x 



6 



(34) 



In order to arrive at the MS scheme, starting from results in the four-dimensional helicity 
scheme given by Eq. (|31~D , we must subtract the counterterm 



2tt 



-p qg (x)+p$(x 



(35) 



After performing this subtraction we have 

,2 1 



Os_ 

2tt 

as 
2tt 



i-Pj. 

2 ' 



~Pqq( 



X 



— Pqq(x) +P, 



x 



k-Pj. 



- (\nx- L)p qq (x) -p { qq \x) 



(36) 



where [. . .] represents all the non-p gg terms in Eq. ([H]). We have used the shorthand 
L = ln(2pj -pj/n 2 ). With the above mass factorization procedure understood, we now 
introduce the notation, 



[dpk{pi,Pj,x)] x(Vfi J ) 



V™ d 5(l-x) + vr + V; 



plus 



(37) 



where we have split the result into an end-point contribution proportional to 5(1 — x) 
regular terms and a term containing 'plus'-distributions. In the q —>■ q case we have, 
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ii,qq 



- + -[--L) + — 



L 2 
2 



1 

2 
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v re9 



■\jplus 

"ii, qq 



1+X 2 

1 - x - In x - (1 + x) (L + 2 ln(l - x)) , 



x 



2L 



l-x) 



+ 4 



ln(l — x) 



x 



The expressions for the other parton splittings, derived in the same manner are, 
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2 1 
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N 

x . \ 21nx 
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x 



2L 
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+ 4 
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(38) 



(39) 



A. 2 Initial-final 



We now consider the case of an initial emitter i with respect to a final state spectator j. 
We first define the dimensionless variables x and u to be, 



X 



The dipole functions are then given by, 
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u if,qq 



if, gq 



g 2 ^ 



ik,j 



f, qg 



ik,j 

if, gg 



XPi ■ Pk V 1 — X + u 

g 2 ^ 

XPi ■ Pk 

g 2 ^ 

XPi ■ Pk 

g 2 ^ 

XPi ■ Pk 



Sik $i 



-1-X 



(l-2x(l-x)), 



X + u 



-2 + 2x(l - x) 



(40) 



(41) 
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with the additional variables needed to account for spin correlations, 



p«j = gV £ 2(1 -x) u(l - u) 

lf ' 9q xpi-pk x Pj-Pk 
V ik f ' j = V ik f ' j , 

if, 99 l /j 91 ' 

gM = Pk_^j_ (42) 
u 1 — u 

The phase-space convolution becomes 

d(f) in) (p h p b -> . . -,Pj,Pk) = J o dx d(j) {n ~ l \p u p b ...,pj) [dpkip^p^x)] , (43) 

where the transformed momenta pi and are defined in terms of the original momenta 
by, 

Pi = xpi, pj = pj +p k + (l- x)pi, (44) 
and pf = pj = 0. The dipole phase space is, 

[dp k (pj; Pi ,x)] = jf^S + (pl) &(x)&(l-x) J— . (45) 

Using the kinematic variables in Eq. fl40|) , the phase space in Eq. fl45|) can be written as 
follows 

[d Pk (Pj;Pi,x)) = ^q^ 2 — du dx Q(u(l - u)) Q(x(l - x)) 

■ (u(l-u)P (1 -x)~ e (46) 

where dfld-3 is an element of solid angle perpendicular to pk and p\. 
We manipulate further as before to yield, 

^ \^,»,*)] = ^%(^)\i-*r*»-^i-»r. (47) 

The integration of the dipoles with this phase space and the mass factorization pro- 
cedure are handled in the same way as the initial-initial case detailed above. The only 
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challenging integral is given by 
% = 



1 i 

du u~ l ~ e (l — u)~ e -j- r 

v ' (1-x + u) 

ir 2 (l-e) (2-x) 



Z_ F (_ 6j _ 2e; 1 - 2e; — - 



(48) 



eT(l - 2e) (1 " v ~~ 7 ~ ~' lc l-x' 

where F(—e, — 2e; 1 — 2e; jj^— ) = l + 0(e 2 (l — x)°). Making the separation of the integrated 
dipoles into different types of contribution yields, (in a notation similar to Eq. (|37D). 
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(49) 



A. 3 Final-initial 

We now consider a final state emitter i, with respect to an initial state spectator j. We 
first define the dimensionless variables x and z to be, 



x = 1 — 



(50) 
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The dipole functions are then given by, 



fl ' qq xpi-pk Kl-x + l-z f ' 

n 2 n 2e 

fi ' gq xpi-p,/ 

n 2 n 2e 
*ik,j y P 



* h 99 xpj -pfoVl — x + 1 — z 1 — X + z J ' 
There is no dipole since the singularities are already accounted for by T>f iiqq . The 

dipole pieces associated with spin correlations are given by, 



-,2,,2e 

1 ' yq xpi -p k Pi- Pk 

fi,gg fi,gq> 



q» = zpt - {1 - z)pl (52) 
The phase-space has the following convolution structure, 

d(j) {n \pa,Vj -> • • • ,Pi,Pk) = / dxd<l>( nr ~ 1 ')(p a ,xpj,-+ ...,Pik) [dp k {puPj,x)} , (53) 



where 

~ d d pk 1 

[dp k (pi;pj,x)] = J^y^ 6+ (Pk) Q(x)Q(l -x) y— 

and the transformed momenta are given by, 



z 



(54) 



Pi = Pi + Pk + (1- x)Pj 

pj = xpj. (55) 
Equation (1S1) can be written out more explicitly as, 



[dpk{pu Pj, x)] = ^ 2P fj\ — ^t-if dz dx e> ( 2; ( 1 ~~ 2 )) ®( x ( l - x )) 



16tt 2 (27T) 1 - 

x (z(l - z))~ e (1 -x)' e , (56) 
where dQj-3 is an element of solid angle perpendicular to the plane defined by pi and pj. 
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The form needed to integrate the dipoles is then, 
■ 2 (IV.,! 3 a s 



9 



[dp k (p ik ]pj,x)} 



x 



(1 - xy x ~ e dzz~\\ - z) 



Pi-Pk' (27T) 1 2e 2ix \2pi -pj 

and the integrations over z can now be performed. The only difficult integral is 
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l fi 
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where F(e, 2e; 1 + e; z) = 1 + 0(e 2 2). With definitions analogous to Eq. fl3T|) we find, 
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(59) 



A. 4 Final-final 



The remaining dipole to consider is one in which we have a final state emitter i with 
respect to a final state spectator j. In this case we define the dimensionless variables y 
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and z, 

The dipole functions are then given by. 



y= , z = ^^. (60) 
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with the auxiliary information, 

n 2 n 2e 9 



//>99 Pi-pkPi-Pk 
Tpkj _ -f^ik,j 

1 fj-!l!l ~ 1 /]■!!■!■ 

# = ^-(i-^pg. (62) 

In terms of the momenta p~j, p~j and p^, this phase-space contribution takes the factorized 
form: 

d(f) {n) (p a ,p + b -> . . . ,Pi,p k ,Pj) = d(p {n - 1 \p a ,p + b^ . . . ,p h pj) [dp k (pi,pj)\ , (63) 
where 

[dpk(Pi,Pj)} = j2n)d=i 5+ (Pk) J(Pk\Pi,Pj) , (64) 
and the Jacobian factor is 

J(Pk;Pi,Pj) = ©(i - *) ©(i - y) (1 7 y)d 3 • (65) 

1 — z 

In terms of the kinematic variables defined earlier, we have 

[dPkfaPj)] = {2P ^ e ^jrk dz dy e(z(l - z)) 8(y(l - y)) 

x (z(l - z)r (1 - y) 1 ^ y~\ (66) 
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where dQ.d-3 is an element of solid angle perpendicular to pij and pu- In this case the 
momenta transform as, 



Pi=Pi+Pk~ 



y 



1 



Pi 



1 ///" " (1 H)'''' 
and the phase-space factorizes to yield the relevant measure, 
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The soft integral is then given by 
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2e 2 T(l - 3e) 

In this case, there is no need to separate the integrated forms of these dipoles and with 
the simplified notation, 



J i d Pk(Pi,Pj,x)} [(T>ff) 



(70) 



we find, 
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Invariant mass distribution in LO and NLO, yu= lOOGeV" 




ni b -b [GeV] 



Figure 3: The bb mass distribution in bb + missing energy events with the renormalization 
and factorization scale \i = 100 GeV. 
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Invariant mass distribution in LO and NLO, // 20GoV 
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ni b -b [GeV] 

Figure 4: The bb mass distribution in bb + missing energy events with the renormalization 
and factorization scale /i = 20 GeV. 
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Figure 5: The scale dependence of the Zbb result at both leading and next-to-leading 
order, with 94 < m b - b < 126 GeV, appropriate for M H = 110 GeV with A = 11 GeV. 
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Invariant mass distribution by pa m or species 
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Figure 6: The contribution of the different parton-parton sub-processes to the cross- 
section, with the renormalization and factorization scale \i = 100 GeV. The lower curves 
represent the leading order results. 
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Figure 7: The significance obtained from our Monte Carlo MCFM, as presented in Table f|, 
compared to the results obtained by the SUSY-Higgs workshop 0. Note that the MCFM 
significance should not be considered realistic because of the lack of detector effects. 
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